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We present a method of disentangling the anisotropies produced by the cliff structures 
in turbulent velocity field and test it in the system of turbulent Rayleigh-Benard (RB) 
convection. It is found that in the RB system the cliff structures in the velocity field are 
generated by thermal plumes. These cliff structures induce asymmetry in the velocity 
increments, which leads us to consider the plus and minus velocity structure functions 
(VSF). The plus velocity increments exclude cliff structures, while the minus ones include 
them. Our results show that the scaling exponents of the plus VSFs are in excellent 
agreement with those predicted for homogeneous and isotropic turbulence (HIT) , whereas 
those of the minus VSFs exhibit significant deviations from HIT expectations in places 
where thermal plumes abound. These results demonstrate that plus and minus VSFs can 
be used to quantitatively study the effect of cliff structures in the velocity field and to 
effectively disentangle the associated anisotropies caused by these structures. 

Key Words: Cliff structures, structure functions, thermal plumes, turbulent thermal 
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1. Introduction 

A paradigm for studying turbulent flow s is the so-called homogeneous and isotropic 
turbulence (HIT) (|Monin fc Yaglomlll975l ). Studying this idealized model allows one to 



focus on the essential physics of small-scale turbulence in the simplest possible case and 
use it as a first step to understand more complicated turbulence problems. However, in 
almost all flow systems existing in nature, anisotropy is always present and unavoidable. 
How to disentangle the effects of anisotropy in the experimentally or numerically mea- 
sured physical quantities has been a major focus in turbulence research in recent years 
and several methods, such as the SO (3) group decomposition, have been put forward to 
separate the isotropic and anisotropic contributions in turbulent flows (Arad et al. 1998, 
1999; Grossmann, VON DER Heydt & Lohse 2001; Biferale et al. 2002). 

Buoyancy- driven turbulent flows occur widely in geophysical and astrophysical sys- 
tems and in numerous engineering applications. In buoyancy-driven thermal turbulence, 
buoyant structures, such as thermal plumes, are predominant coherent structures that 
transport heat and drive the flow (Shang et al. 2003; Xia, Sun & Zhou 2003). Earlier 
visualization experiments (Moses, Zocchi & Libchaber 1993; Xi, Lam & Xia 2004) have 
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shown that these structures consist of a cap with sharp temperature gradient and a 
stem that is relatively diffusive and hence would generat e cliff-ramp-like structures in 
temperature time series when passing a thermal probe ( Belmonte fe Libchaber 19961 : 
Zhou fe Xiall2002h . It is well known that the so-called cliff-ramp structur es would induc e 
strong anisotropic effects. This has been widely studied in passive scalars (| WarhaftfeOOOl ) . 
but the study on the effects of buoyancy-induced anisotropy is very limited. Here, we use 
turbulent Rayleigh-Benard (RB) convection, a fluid layer heated from below and cooled 
on the top, as an example to study the anisotropic effects induced by buoyancy. In the 
past few decades, turbulent RB convection has become a model system for studying the 
phenomena and the generic physics associated with turbulent flows driven by buoyancy 
(Ahlers, Grossmann & Lohse 2009; Lohse & Xia 2010). 

In the field of turbulence, the velocity structure functions (VSF), namely 



Sp (r) = (\S r vn, 



(1.1) 



are usually used to characterize the turbulent kinetic energy cascades, and hence they 
are of prime importance and have been the central focus in the study of fluid turbulence 
(Sreenivasan & Antonia 1997; Ishihara, Gotoh & Kaneda 2009). Here, 5 r v — v{x + r) — 
v(x) is defined as the velocity increment over a separation r and (• • •) denotes an ensemble 
average. In part icular , the situa t ion fo r thermally-driven turbulence is more complicated. 



Bolgianol ( 19591) and lObukhovl ( 19591 ). denoted hereafter as B059 for short, have long 



argued that within the so-called inertial range and above a certain buoyant scale, i.e. 
the Bolgiano length scale £b, buoyant forces drive the cascade processes and scale as 
S p (r) ~ r 3p / 5 . However, experimentally whether the BQ59 scaling exis ts in turbulent 
RB system remains unsettled (for a recent review, see lLohse fe XialboiOl ). The maturing 
of spatial- velocity- field measuring techniques, e.g. the particle image velocimetry (PIV), 
has provided a great impetus for our understanding of such issues (Sun, Zhou & Xia 2006; 
Zhou, Sun fc Xia , 2008; Kunnen et al. 2008; Zhou & Xia 2010). In particular, the work of 
Sun et al. (2006) have shown that, at the center of the convection cell, the velocity field 
exhibits the same scaling behavior that one would find in HIT, whereas, near the cell 
sidewall, a "new scaling" was found. While results at both places imply no B059, the 
observed scaling behavior near sidewall remains unexplained by any existing theoretical 
models. We note that in a closed RB cell, the spatial distribution of the buoyancy-driven 
thermal plumes is hi ghly inhomogeneou s , they abound nea r the s i dewall but are scarcely 
found at the center (|Qiu fc Tonel EoOll IShang et aZ.I I2003L booi IXi et all 12004 ) . These 
give us a clue that plumes may be responsible for the different scaling behavior observed 
at different places of the system and suggest that buoyancy plays an important role 
in the cascade process, but how this comes about is still missing. The objective of the 
present study is to address such a question, i.e., how buoyant forces influence the cascade 
properties in turbulent RB system? 



2. Experimental setup and parameters 

The convection cell has been described in detail elsewhere ( Zhou fc XiaboiOl :! IZhou et al. 



20111 ) and here we give only its main features. It is a vertical cylinder of height H = 50 cm 



and inner diameter D = 50 cm and thus of unity aspect ratio. The top and bottom plates 
are made of 1.5 cm thick pure copper with nickel-plated fluid-contact surface and the 
sidewall is made of a plexiglas tube of 5 mm in wall thickness. Deionized and degassed 
water was used as the convecting fluid. A square-shaped jacket made of flat plexiglas 
plates and filled with water is fitted to the outside of the sidewall, which greatly reduced 
the distortion effect to the PIV images caused by the curvature of the cylindrical sidewall. 
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Figure 1. Left: A snapshot of the instantaneous vertical velocity field w(x,z) near the cell 
sidewall for Ra — 4.0 x 10 10 . Positive is defined for the upward motion and color coding is in 
cm/s. Right: Vertical slices of the vertical velocity component w(x,z) through the left figure. 
The red curves mark the extracted cliff structures, corresponding to the regions of plume fronts. 



Two series of measurements of the spatial velocity field were carried out using the PIV 
technique. In the first series, the measuring positions were fixed near the cell sidewall and 
the experiments covered the range 5.9 x 10 9 < Ra < 1.1 x 10 11 of the Rayleigh number 
Ra = agATH 3 /vk, with g being the gravitational acceleration, AT the temperature 
difference across the fluid layer, and a, v and k being, respectively, the thermal expansion 
coefficient, the kinematic viscosity, and the thermal diffusivity of water, whereas in the 
second series, the measurements were made from near the cell sidewall to the cell center 
at fixed Rayleigh number Ra = 4.0 x 10 10 . For both series, the cell was tilted by a 
small angle of about 0.5° so that both series were made within the vertical plane of the 
large-scale circulation and at midheight of the cell. During the experiment the entire cell 
was wrapped by several layers of Styrofoam and the mean temperature of the convecting 
fluid was kept at 29°, corresponding to a Prandtl number Pr = v/k = 5.5. The details 
of the PIV measurement could be found elsewhere (Xia et al. 2003; Sun, Xia & Tong 
2005), here we give only its main features. Hollow glass spheres of 10 /im in diameter 
were chosen as seed particles and the thickness of the laser lightsheet was ~ 0.5 mm. The 
spatial resolution of the measured velocity field is 0.59 mm, which is much smaller than 
the lower end of the inertial range and hence is sufficient to reveal the scaling properties 
in the inertial range. In each measurement, the measuring region has an area of 4.7 x 3.7 
cm 2 (see the left panel of figure [T]), corresponding to 79 x 63 velocity vectors, and the 
experiment lasted 3 hours in which a total of 25000 two-dimensional vector maps were 
acquired with a sampling rate ^2.3 Hz. Since buoyant forces are exerted on the fluid in 
the vertical direction, we focus our attention mainly on the longitudinal vertical velocity 
increments S r w = w(x, z + r) — w(x, z), where w(x, z) is the vertical velocity component 
obtained at position (x, z). To acquire the accurate statistics, both temporal average and 
spatial average of moments of velocity increments within an area of 1.8 x 3.7 cm 2 , i.e. 
31 x 63 velocity vectors, were used when calc ulating VSFs, as the flow is approx imately 
locally homogeneous in turbulent RB system ( Sun et al. 20061 : Zhou et al. 20081 ). 
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Figure 2. (a) PDFs of the normalized increments of vertical velocity components at x/D = 0.04. 
From the inner to the outer PDF, r/r\ — 16.6 (dark-green up-triangles), 33.3 (pink squares), 
57.1 (blue down-triangles), and 107 (red circles), all within the inertial range, (b) Skewness of 
8 r w as a function of r jr\ for different measuring positions x/D from near the cell sidewall to the 
cell center. The data were obtained at Ra = 4.0 x 10 10 . (c) Skewness of 5 r w as a function of 
r/n for five different Ra obtained near the cell sidewall (x/D — 0.04). 



3. Results and discussions 

3.1. Cliff structures in the vertical velocity field 

The left panel of figure 1 shows a typical snapshot of the vertical velocity field w(x, z) near 
the cell sidewall, where x is the horizontal distance from the wall. The arrow in the figure 
indicates the region with large vertical velocity, which is typically caused by a group of hot 
plumes passing by. The plumes are generally believed to be detached thermal boundary 
layers by buoyant forces. As introduced in § 1, thermal plumes usually generate cliff 
structures in the temperature field. Here, the snapshot of w(x, z) and the associated slices 
through w(x, z) in the vertical direction (the right panel of figure [1} further show that 
in addition to temperature, plume fronts can also generate cliff structures in the vertical 
velocity field. The formation of cliff structures in the velocity field may be understood 
as follows: Under the action of buoyant forces thermal plumes possess a higher speed 
in the vertical direction in comparison to the background fluids, which would deform 
the plumes such that they are compressed in the vertical direction and stretched in the 
horizontal directions. This deformation shortens the distance between the (relatively) 
high speed plume front and the low-speed downstream fluids, therefore resulting a steep 
velocity gradient, i.e. a cliff structure. 

One can expect that the persistence of such cliff structures would violate the local 
isotropy of turbulent flows. Isotropy is the central hypothesis for most theories and models 
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of small-scale turbulence. Although a number of investigations in the atmosphere and 
in laboratory flows have sho wn that the skewness of velocity derivative is nontrivial 



( Sreenivasan fc Antonialll997l ). implying local anisotropy at small scales, the presence of 
cliffs would enhance the degree of this small-scale anisotropy. This can be reflected by the 
probability density functions (PDF) of velocity increments over different scales. Figure 
[2{a) shows PDFs of 6 r w normalized by the standard deviation of w, a w — J ((w — (w)) 2 ), 
over several different length scales within the inertial range. The data were obtained near 
the cell sidewall (x/D = 0.04) at Ra = 4.0 x 10 10 . The asymmetry of the distributions 
with long left tails is clear, especially for the (relatively) larger scales, which signifies the 
anisotropy at these scales. To quantify this anisotropy, we plot in figure[2j6) the skewness 
of S r w as a function of r/?70 for different measuring positions ranging from x/D = 0.04 
(near the sidewall) to x/D — 0.50 (at cell center). It is seen that all skewness are negative, 
but their magnitudes decrease continuously when moving away from the wall and appear 
to remain invariant for x/D > 0.4. If the nonvanishing skewness observed here is indeed 
induced by cliff structures of thermal plumes, the behaviors of the skewness may then 
be understood from the inhomogeneous spatial distributions of thermal plumes in the 
convection cell, i.e. , plumes abound near th e sidewall but a re scarce in the central region 



( Qiu fc TonjbOOlUShang et aLll2003L [20041: IXi et al\\2004) 



Figure^c) shows the skewness of S r w as a function of r/rj for five different Ra obtained 
near the cell sidewall (x/D = 0.04). Again, all skewnesses are found to be negative. An- 
other noticeable feature is that the magnitudes of the skewness decrease with increasing 
Ra. This may be understood as the cliff structures are being smoothed out by the in- 
creased turbulent fluctuations. As the cliff structures of the velocity field are generated by 
buoyancy, this result suggests that temperature becomes more passive when the convec- 
tive flow becomes more turbulent. Note that previous experiment al study has shown that 



temperature may becom e progressively passive for Ra > 10 ijBelmonte fc Libchaber 



19961 : IZhou fc Xiall2002h . which is qualitatively consistent with the picture obtained in 



the present study. 

3.2. Plus and minus velocity structure functions 

The fact that plumes can induce anisotropy suggests that the cascades of velocity field 
would possess different dynamics in different directions, i.e., the positive and nega- 
tive velocity increments may have different scaling behaviors. To study this quantita- 
tively, we examine the statistical properties of the plus and minus longitudinal VSFs 
( Vainshtein fc Sreenivasan 1994 ISreenivasan et aLlll996h . defined as 



S±(r) = {[(\5 r w\±5 rW )/2F). (3.1) 

From this definition, it is clear that cliff structures such as those shown in the right 
panel of figure [T]) will be excluded from «S+(r) but included in S~(r). This enables one 
to study separately the contributions of positive and negative velocity increments to the 
corresponding structure functions. A similar analysis has been performed previously in 
a turbulent channel flow and it was found that the plus VSFs are less affected by the 
presence of the wall ([Onorato fc Iusdl200lh . Figure a) plots in log-log scale the com- 



pensated third-order VSFs S^(r)/r (triangles), Ss(r)/r (circles), and S^(r)/r (squares) 
vs r/rj near the cell sidewall, which exhibit slightly different scaling ranges. The compen- 
sated plot shows a flat range for S^(r)/r, i.e. S3 (r) ~ r, suggesting that S3" (r) possesses 

f The Kolmogorov length scale n is estimated from n = (^ 3 /e) 1//4 , where e = o%/L is the 
energy dissipation rate per unit mass and L is the largest length scale of the turbulence. It is 
found that n changes from 0.27 mm near the sidewall to 0.26 mm at cell center for Ra = 4.0 x 10 10 . 
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Figure 3. (a) Compensated third-order VSFs S$(r)/r, Ss(r)/r, and S%(r)/r measured near 
the cell sidewall. The arrow marks r = 1.5^3 for reference, (b) Comparison of VSF exponents 
£p , Cpj an d Cp~ with various model predictions. The data were measured at x/D = 0.04 for 
Ra = 4.0 x 10 10 . (c) The third-order VSF exponents £+, £ 3 , and C 3 " 

as a function of the 



normalized distance from the cell sidewall x/D for Ra = 
Qi, and ^ obtained near the cell sidewall (x/D = 0.04). 



4.0 x 10 . (d) i?a-dependence of C3 



the same scaling behavior as that for HIT, whereas both Ss(r) and S^ir) exhibit much 
steeper scalings. To reveal this more clearly, we show in figure the measured scal- 
ing exponents £+, ( p , and CjT of the VSFs of o rders p = 1 to 8. Wh en comparing 
with the predictions of the hierarchy models of She &: Leveauel (|l994h (SL94) for HIT, 
we find excellent agreement. These results suggest that the scaling behaviors of S+(r) 
are consistent with what one would expect for HIT. Moreover, because cliff structures 
contribute mainly to the negative velocity increments, £~ should exhibit some deviation 
from the K41-type scaling. This is indeed observed. Figure^ b) shows that £~ and ( p are 
both much larger than the predictions of SL94. Therefore, the plus and minus structures 
functions are effective means to study the effects of cliff structures in the velocity field 
and can be used to effectively disentangle the associated anisotropies caused by these 
structures. These cliff structures are produced by thermal plumes in the present case, 
but in general they can be produced by coherent structures in other types of flows. 

Let's move on now to the location- and i?a-dependencies of the VSF scaling exponents. 
We focus mainly on the third-order because £3 = 1 is an exact result for HIT. Figure 
[3fc) shows the exponents C3" j (3, and C3 as a function of the measuring position x/D. 
It is seen that both (3 and C3 are much larger than the K41-value of 1 near the cell 
sidewall and drop from the wall. For x/D > 0.3 £3 and C3" are both essentially 1. On 
the other hand, apart from some data scatter, (3 assumes essentially the K41 value for 
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Figure 4. Comparison of ESS VSF exponents Cp3i £~ 3 , and the SL94 scaling exponents. 



all x/D. The decreases of £3 and £3" correspond to the reduced anisotropy associated to 
the cliff structures, i.e., the number of plumes decreases as the measuring position moves 
away from the cell sidewall, confirming quantitatively the results shown in figure 
For x/D > 0.3, i.e. the cell's central region, thermal plumes are scarce and hence it is 
not surprising to obtain approximately the K41 scaling for all three VSFs. 

Figure OUJd) shows the measured scaling exponents versus Ra. Again, £3" is seen to 
remain nearly constant around the value of 1, but (3 and £3", despite certain scatter, 
show an overall decreasing trend with increasing Ra. Here, the decrease of £3" implies 
that the influence of buoyancy on the cascade processes becomes weaker when the flow 
becomes more turbulent, which could also be reflected by the behavior of the skewness 
of 5 r w, whose magnitude is found to increase with decreasing Ra as shown in figure^c). 

It should be noted that near the sidewall the result ( 3 > £3", shown in figurc[31 does not 
imply that the positive velocity increments are more intermittent than the negative ones. 
To illustrate thi s, we examine the scaling behavior of VSFs via the extended self-similarity 
(ESS) method (jBenzi et aZ.lll993h . i.e., S p (r) is plotted against S3 (r), instead of r, in a 
log- log scale. Figure [4] shows the measured ESS (relative) scaling exponents £+ 3 and (~ 3 
for the plus and minus VSFs, respectively. One sees that (~ 3 are slightly smaller than 
(p 3 . This result suggests that it is the minus velocity increments that possesses a higher 
degree of intermittency, which should be attributed to the persistence of cliff structures. 
A detailed analysis of such cliff structures would therefore be helpful for understanding 
the present results. To do this, we use a criterion to identify cliff structures i n the vertical 



slices of w(x,z) which i s similar to those used for passive and active scalars (IMoisv et 
l200lUZhou fc XiakOOdh a cliff is ident ified when —dw/dz > ct w /Ib, where £3 is based on 



the global quantities ( Sun et a"Zll2006l) . When such a cliff is found, we define its position 
zq as the point maximizing \dw/dz\, and its width Ac as the separation in space between 
the two extrema of w surrounding z . Applying this procedure, over 50 000 cliffs were 
identified near the sidewall (x/D = 0.04) for each Ra. Three examples of the extracted 
cliff structures are shown as red curves in the right panel of figure [T] Figure [S] shows the 
mean cliffs' width (Ac), normalized by Ib- One sees that (Ac) = (1.5 ± 0.2)£b is nearly 
independent of Ra. We further note that the cliff structures mainly occur at the scales 
near the upper end of the VSF scaling range [see figure |3]J a )} - This could also be reflected 
qualitatively from figure [2j a), which shows that the distributions of S r w over large scales 
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Figure 5. Mean cliffs' width (Ac), normalized by the Bolgiano length scale £b, as a function of 
Ra obtained near the cell sidewall (x/D — 0.04). The error bars mark the standard deviation 
of Ac- 



seem to be more asymmetric than those over relatively small scales, and quantitatively 
from figure from which one sees that within the inertial range (16 < r/rj < 110) 

the magnitude of the skewness of S r w increases with the scale r for the near wall data 
(x/D = 0.04, dark-green up-triangles) . Cliff structures contain large velocity variations 
and would enhance the magnitude of velocity increments over the scales around their 
mean width (Ac), i-e. the scale near the upper end of the VSF scaling range. This can 
also be seen from figure a) that and S~ differ the most around 1.5£b, signifying 
that scale is representative of the typical size of a cliff structure. The cliffs' influences 
to the scales near the lower end of the scaling range, however, are much weaker [see 
figure [3j a)]. Therefore, under the influences of cliff structures, the value of Sp(r) near 
the upper end of the VSF scaling range would increase, while those near the lower end 
do not change significant. As a result, the measured scaling exponents of S~(r) would 
increase. With decreasing buoyancy effects, the scaling exponents are also expected to 
drop, which is indeed observed in figure EJd). 



4. Conclusion 

To summarize, we have demonstrated that as a manifestation of buoyancy effects ther- 
mal plumes can generate cliff structures in the vertical velocity field, which would in 
turn generate asymmetry in the velocity increments. We further show that such effects 
can be quantified by examining respectively the plus and minus velocity increments. For 
the plus increments, which largely exclude the cliff structures and hence removing the 
buoyancy effects, the scaling of their moments are the same as those expected of homo- 
geneous and isotropic turbulence (HIT). For the minus increments, the scaling is found 
to be much steeper than that for HIT, which is caused by the presence of cliff structures. 
Such effects of buoyant forces are found to vanish gradually when moving away from the 
cell sidewall, owing to the inhomogeneous distributions of thermal plumes in a closed 
convection cell. It is also shown that, as a result of cliff structures being smoothed out by 
the increased turbulent fluctuations, the effect of buoyancy on the velocity field decreases 
with increasing Ra. Despite its simpleness, the analysis presented here provides a useful 
method to quantitatively study the effect of buoyancy and disentangle its contribution 
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in the velocity field in buoyancy-driven turbulence. More generally, it may be used to 
quantify the effect of anisotropy in other complex phenomena when sharp fronts exist in 
the related random fields. 
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